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SCALING LIMITS OF SOLUTIONS OF LINEAR EVOLUTION 
EQUATIONS WITH RANDOM INITIAL CONDITIONS 


MILOSZ KRUPSKI 

Uniwersytet Wroclawski, Instytut Matematyczny 
pi. Grunwaldzki 2/4, 50-384 Wroclaw 

Abstract. We consider a linear equation dtu = On. where £ is a generator 
of a semigroup of linear operators on a certain Hilbert space related to an 
initial condition u((Y) being a generalised stationary random field on R d . We 
show the existence and uniqueness of generalised solutions to such initial value 
problems. Then we investigate their scaling limits. 


1. Introduction 

Partial differential equations with random initial data appear in a wide variety 
of topics, e.g in the study of the Large Scale Structure of the Universe or growing 
interfaces in deposition of chemical substances [AMS94I , iJWOll , |Woy98| . 

In particular, such equations have been studied in a number of papers which were 
primarily concerned with calculating scaling limits of solutions LW981 . AL02] . 
IRMAAOll . AL01I . However, they all lack a formal definition of solutions and 
instead depend only on classical formulas for non-random equations, extended to 
accommodate random fields. Furthermore, the absence of a proper definition pre¬ 
vented the authors from exploring conditions for the existence and uniqueness of 
solutions. 

Moreover, in the context of scaling limits of solutions, it is a welcome result to 
show that the obtained limit (provided it exists) is also a solution to some problem 
related to the original one, perhaps only in a less regular sense. Looking at previous 
results, e.g. ILW981 , |AL021 . one may indeed see this being the matter; however in 
most cases, the concept of a generalised random field is necessary to properly define 
initial conditions of those derived problems. 

In this paper we put forward a comprehensive approach to study general linear 
evolution equations supplemented with random initial conditions, using tools pro¬ 
vided by the theory of generalised random fields. We consider the following abstract 
initial value problem 

j d t u = Cu on [0, oo) x R d , 

y u(0) = 770 on R d , 

where 770 is a broad-sense stationary generalised random field and £ is a linear 
operator generating a Co-semigroup of linear operators on a certain Hilbert space 
corresponding to 770 . 

In Section 2 we prepare tools necessary to develop our theory. Section 3 contains 
main results, including the definition of solutions and a proof that, under reason¬ 
able assumptions, a kind of a “semigroup solution” is in fact unique. Section 4 is 
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concerned with calculating scaling limits of solutions in an abstract setting. Finally, 
Section 5 contains examples illustrating the results. 

The scope of this work is limited to linear problems; however, the ideas it contains 
seem to be applicable to non-linear problems as well. Such considerations will be 
the subject of our forthcoming papers. 

2. Preliminaries 

2.1. Basic notation. We denote the Borel sigma-algebra on by Bor(K. d ) and 
the Lebesgue measure by either A or dx. 

We use the Fourier transform defined as 

J 7 {f}(0= [ e*- x f{x)dx. 

J ! d 

Given a measure space (X, 0,/z), by L 2 (X, Q,p) we denote the space of all 0- 
measurable real functions such that the integral 

[ \f?dn= [ \f(x)\ 2 n(dx) 

Jx Jx 

is finite. Usually we will shorten the notation to L 2 (/z) = L 2 (X, 0, p). 

Let us fix a probability space (U, E, P) and denote L 2 (U) = L 2 (U, E, P). This is 
the set of all random variables with bounded variations. It is a Hilbert space with 
the standard inner product EXY defined for all X, Y £ P 2 (U). 

We also use the notation L 2 (X, 0,/x; C), P 2 (U;C) to describe analogous spaces 
of complex-valued functions. 

We write 

(Xi,..., X fc ) = (Yi,..., Y*,), where X i: Y { £ P 2 (U) and k £ N 
if both random vectors have the same probability distributions and we say that 
w-lim(X 1 „,..., Xfc „) = (Yi,..., Yk) 

n—too 

if the corresponding probability distributions converge weakly as measures on the 
product space. 

For a linear operator A, by Dom(A) we denote its domain. Let us also recall a 
definition of an admissible subspace following IPaz831 Ch. 4, Def. 5.3]. 

Definition 2.1. Let Y be a linear subspace of a Banach space (X, || • ||.y) such that 
there exists a norm || • ||y in which Y is a Banach space itself, and ||y||jc < Cj|x/||y 
for all y £ F. Suppose A is the generator of a Go-semigroup T(t) on X. A subspace 
Y of X is called Al-admissible if it is an invariant subspace of T(t), t > 0, and the 
restriction of T(t) to Y is a Co-semigroup on Y (i.e. it is strongly continuous in the 
norm || • ||y). 

2.2. Generalised random fields. Now we recall the basics of the theory of gener¬ 
alised random fields. For a more detailed description we refer the reader to [ QV64L 
IB.oz981 and llJrb58j . 

First we introduce the space of test functions D = C^° (R d ) consisting of smooth 
functions with compact supports. 

Definition 2.2. A sequence (p n is convergent to <p in D if 

(1) there exists a compact K C R d such that supp ip U U^Lo SU PP <Pn C K 

(2) derivatives d a tp n converge uniformly to d a tp for every multi-index a. 

Definition 2.3. A generalised random field is a continuous linear operator y : D —> 
L 2 (U). 
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With a slight abuse we use the following notation 
vivi, ■■■,<Pk) = (ri(.<Pi), ■ ■ ■ > v(<Pk))- 

For two generalised random fields we write 77 = £ and say that 77 and £ are equiv¬ 
alent in finite-dimensional distributions if for every k £ N and every vector of test 
functions (ipi,... ipk) £ D k we have 

ri(tpi,...,ip k ) =£(<Pi 

Similarly we say r] n converge weakly in finite-dimensional distributions to 77 if for 
every vector of test functions (ipi,... ifk) £ D k 

w-lim 77 „ (y>i, ...,ip k )= vivi, ■■■, Vk), 

n—> 00 

and write w-linin^oo r) n = 77 . 

2.3. Random measures and integrals. 

Definition 2.4. Let a be a Borel measure on W 1 and II C Bor(R d ) be a family 
of all sets A such that cr(A) is finite. An orthogonal random measure Z with a 
reference measure cr is a function Z : II —> L 2 { 12; C) such that 

(1) EZ(A) = 0, 

(2) EZ(A 1 )Z(A 2 ) = a(A 1 nA 2 ), 

(3) ^(Un=i A n ) = ( as a limit in L 2 (fl;C)) for every pairwise 

disjoint collection of sets {A n } C II such that (J A n £ II. 

Notice that, despite the name, in general Z is not a measure since II is not 
necessarily a sigma-algebra. It may even be the case that Z cannot be extended 
to a proper measure. It is however possible to define an integral with respect to Z 
(see e.g. |Kry02[ Ch. 2, Sec. 3] or IGS691 Ch. 5, Sec. 3] for more details). Here, let 
us only briefly recall that the properties of an orthogonal random measure allow us 
to define a unique linear operator Iz : L 2 (R d , Bor(R d ), cr; C) —> L 2 (H; C), such that 

(!) Me Cil Ai ) = E CiZ(Ai) a.e. for all n £ N, Ai £ n and c, £ C, 

i<n i<n 

and 

E\Iz{f) r= / \f\ 2 da for all / £ L 2 (R d , Bor(R d ), cr; C). 

J R d 

Remark 2.5. In order to define an othogonal random measure it suffices to con¬ 
sider a family of sets n 0 such that Ai,A 2 £ n 0 implies A\ fl A 2 £ n 0 (a tt- 
system) and there exist A\ C A 2 C ... £ n 0 such that ( A n = (see |Kry02[ 
Ch. 2, Sec. 3, Thm. 19]). 

We say that Z n converges weakly in finite-dimensional distributions to Z if for 
every vector of functions (/ 1 ,..., fk) £ L 2 (R d , Bor(R d ), cr; C) k we have 

w-lim (IzMi), ■ ■ = (Wil ■ ■ -,Wk)), 

n—too 

denoting w-lim^oo Z n = Z. 

Proposition 2 . 6 . Suppose we have 

lim E\Z n (A) — Z(A )| 2 = 0 for every A £ II. 

n—too 


Then w-lim. 
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Proof. Because of © and the continuity of the operators Iz n and Iz, we have 
lim E\IzM) - !zU)\ 2 = 0 for every / £ L 2 (K d , Bor(R d ), cr; C). 

n—> oo 

Denote 

k n =E\I Zn (f)-Iz(f)\ 2 . 

Then by the Chebyshev inequality we have 
P{\Iz n {f) ~ Iz{f)\ >e)< 

which shows that for every / £ L 2 (R d , Bor(ffi. d ), cr; C), Iz n (f) converges in proba¬ 
bility to Iz(f) and therefore also 

w-lim I Zn (f) = Iz(f)- 

n—> oo 

Then by linearity of the operators Iz n and Iz and the Cramer-Wold theorem 
(see [Bil951 Thm. 29.4]) we obtain w-limn-^ Z n = Z. □ 

From now on we shall write 

Iz(f)= f f{x)Z{dx). 

J R d 


2.4. Broad-sense stationary random fields and H space. For tp £ D and a 

parameter h £ R d we denote the translation operator by 

r h ip{x) = <p(x + h). 

Definition 2.7. A generalised random field rj is broad-sense stationary if for every 

ip, i/j £ D 

Ei= Erj(ip), 

Er]{ ( P) r ](' l l ; ) = Erj(Thp)ij{Thijj) for every h £ K d . 

In this work we limit ourselves to broad-sense stationary random fields with zero 
mean, i.e. Eri(p) = 0 for every <p £ D. 

We recall two basic theorems from the theory of broad-sense stationary random 
fields, the proofs of which may be found in IGV64] or |GS69] . A Borel measure a 
such that 

J R d (1 + \x\ 2 Y 

is finite for some p > 0 is called tempered. 

Theorem 2.8. Let p be a broad-sense stationary random field. There exists a 
tempered measure a, such that 

Er){<p)r)(ip) = [ E{fi}E{p}dc r. 

The measure a is called the spectral measure of the random field p. 

Theorem 2.9. Let p be a broad-sense stationary random field with the spectral 
measure cr. There exists an orthogonal random measure Z with the reference mea¬ 
sure cr such that for every ip £ D 

(2) p(p)= f T{p\dZ. 

J R d 
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Conversely, given an orthogonal random measure Z , by the means of formula (121), 
we may construct a generalised broad-sense stationary random field. 

Using the representation from Theorem 12.91 we may also consider the following 
continuous extension of a broad-sense stationary random field 77 to a space larger 
than D. 

Given the spectral measure a associated with 77 , we consider the space D as a 
pre-Hilbert space with the inner product 

(3) ( f,g)=[ E{f}X{g}da + [ fgdx. 

J R d J R d 

We define the space H (a) as the completion of D with respect to ([3]) . Then H(a) is 
a Hilbert space with the norm = (/,/)■ Furthermore, D is a dense linear 

subspace of H(a) and thus H{a) is a dense linear subspace of i 2 (R d , Bor(R d ), da;). 
We may now define an extension 77 : H(a) —> L 2 (Q) by the formula 

v(f) = [ Hf) dZ - 

J R d 

Since for every / £ H(a) we have 


EW)) 2 = e 



F{f}dZ 


[ \E{f}\ 2 da< \\f\\ 2 H{a) < 00 , 
JR d 


this extension is well-defined and continuous. 


Remark 2.10. H(a) is a separable space. 

2.5. Transformations of stationary random fields. Consider a broad-sense 
stationary random field 77 with spectral measure a , extended to the space H(cr) as 
described in the previous section. Let A : L 2 (dx) —> L 2 {dx) be a continuous linear 
operator. Its standard adjoint operator A* is determined by the formula 

f Afgdx= f fA*gdx for all /, g £ L 2 (dx) 

Js. d J R d 

Definition 2.11. Let A : L 2 (dx) —> L 2 (dx) and assume H(a) is an invariant 
subspace of A*. A field Ai 7 : H(a) —> L 2 (VL) defined as 

Ar](<p) = r](A*<p) for every ip £ H(a) 

is called a transformed broad-sense stationary random field. 

Remark 2.12. Let us emphasise that the operator A*, in general, is not the adjoint 
of A on H(a). 

Definition 2.13. Suppose rj f are broad-sense stationary random fields with spectral 
measures at and consider a Banach space X such that X C H(at) for every t> 0. 
A time-dependent family At r]t of transformed broad-sense stationary random fields 
is called locally integrable with respect to X if for every if) £ C c ([0, 00 ), X) 

( 1 ) t !->■ At ? 7 t(^(t)) is a measurable function, 

(2) f 0 °° E(At7]t(ifi(t))) 2 dt < 00 . 

An immediate consequence of this definition is that the mapping 1 1 —> Atiit{il)(t)) is 
Bochner-integrable (see 1DU771 Ch. 2, Sec. 2, Thm. 2]). 

2.6. Regular stationary random fields. Having defined generalised random 
fields, we should also mention their non-generalised, regular counterparts. We re¬ 
strict ourselves to the one-dimensional, broad-sense stationary case for simplicity 
of the exposition. 
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Definition 2.14. A continuous mapping 77 : M — > L 2 (Q) is called broad-sense 
stationary if for every x, y € R 

Ev(x)v(y) = R{\x-y\), 

for some R : R —> R. 


ft is easy to see that R must be continuous and positive-definite, thus the 
Bochner theorem asserts that there exists a finite measure cr on R, such that 
R(x) = f R e* x f cr(d£). In full analogy to Theorem 12.91 we also have the representa¬ 
tion r](x) = f R e lx Z Z{d^), where Z is an orthogonal random measure with reference 
measure a. 

Since 

/ \T{tp}\ 2 da < cr(R)(A(supp ip) sup \tp\ ) 2 < 00 , 

Jr 

we may define the corresponding generalised broad-sense stationary random field 
77 0 by using the formula 


V°(<P) = [ dZ. 

Jr 

Conversely, suppose we have a generalised broad-sense stationary random field rj °, 
with its associated orthogonal random measure Z and a finite spectral measure cr. 
Then E(f R e^Z(dO) 2 = cr(R) and we may define a regular broad-sense stationary 
random field 77 by the formula 



Example 2.15. Consider the standard (two-sided) complex-valued Brownian mo¬ 
tion B t on R and define Wo([0, t}) = B t . We may then extend Wo to an orthogonal 
random measure W on R, defined on the family of bounded sets (i.e. of finite 
Lebesgue measure, cf. Remark l2.5l) . Indeed, we may also check that 

EW([ 0, f])W([0,s]) = EB t ~B~ s = min{s, t} = A([0, f] n [0, a]), 

which proves that A is in fact the reference measure of W. 

Define a generalised stationary random field 77 ( 175 ) = j R R{(p}W(dx). Then, 

Er]{p)r](ip) = f J-{<p} J-{ij)} dx = (2n) d f dx 

Jr Jr 

(the latter equality follows from the Parseval identity). The field 77 is called the 
(Gaussian) white-noise, ft follows immediately that H{dx) = L 2 {dx). Since A is 
not a finite measure, one may show that this random field cannot be represented 
as a regular random field. 


3. Solutions to Cauchy problem 

The theory recalled in the previous section allows us to define solutions to the 
initial value problem 

f d t u = Cu on [0, 00 ) x R d , 

( 4 ) , 

[ it(0) = ? 7 o on R d 

and to study their properties. We assume the following 

(HW: 770 is a generalised broad-sense stationary random field; 

(0)3): £ is a generator of a Co-semigroup {e tC }t>o of linear operators on 
L 2 {dx) such that D C Dom(£*) and H{<jq) is £*-admissible (see Defini¬ 
tion 12 . 11 ) . 
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Remark 3.1. The operator C* is also a generator of a Co-semigroup on L 2 (dx) and 
(e tc )* = e tc (see IPaz831 Ch. 1, Cor. 10.6]). 

Definition 3.2. A time-dependent family u(t) of transformed broad-sense station¬ 
ary random fields, which is locally integrable with respect to H(<jf) (see Defini¬ 
tion [5T3]) is called a solution to problem J]} if for every vector (ip\,... ,ipk) £ 
C~([ 0 ,oo) x R d ) k we have 
/»00 

( 5 ) / u(t)((d t dt = —770(^1 ( 0 ),..., ^fc(O)). 

Jo 

Remark 3.3. Notice that the expression under the integral in (O is well-defined 
since for every ip(t) £ C£°([0, 00 ) x R d ) a direct calculation shows 

(d t + C*)ip(t) £ C c ([0,oo),H(a 0 )). 

In the following theorems we extend the meaning of classical semigroup solutions 
to cases where the initial conditions are broad-sense stationary random fields. 

Theorem 3.4 (Existence of solutions). Assume (|3Ji) and (JDo). The family of 
transformed broad-sense stationary random fields e tc r/o (see Definition \2.11\) is a 
solution to problem ©■ 

Proof. First, we show that the function t 1 —> e tc ip(t) is continuous for every ip £ 
C c ([0,oo),H(ao)). Recall that there exist constants M > 1 and to > 0 such that 
||e* £ \\h(<t 0 ) < Me tul (see [Paz831 Ch. 1, Thm. 2.2]). Denote by T the upper bound 
of the support of ip and let t n —> t. Then 

II e tc *ip(t) - e tnC *ip(t n )\\ H[cro ) 

< ||(e‘ £ * -e*» £ >(*)llff ( «r o) + \\^ c \m-^(tn))\\H {ao ) 

< ll(e t£ * - e tnC *)ip(t)\\ H (ao) + Me Tuj \\ip(t) - ip(t n )\\ H{ao) . 

It follows from the continuity of ip that lim^oo ||^>(f) — ip(t n )\\H{a 0 ) = 0. Since the 
semigroup is strongly continuous, the function t 1 —> e tc p is continuous for every 
ip £ H(ao) (see [Paz83[ Ch. 1, Cor. 2.3]) and therefore 

lim \\(e tC * - e tnC *)ip(t) || ff(ffo) = 0 . 

n—¥ 00 

Because 770 is a continuous operator, the function e tc ri 0 (ip(t)) = rjo(e tC ip(t)) is also 
continuous and hence measurable. Now, 

/•OO /»oo 

/ E(e tc r] 0 (ip(t))) 2 dt = / E(r]o(e tC *ip(t))) 2 dt 

Jo Jo 

= [ [ \E{e tc *ip(t)}\ 2 da 0 dt <T ■ sup \\e tC * ip(t)\\ 2 H(ao) < 00 . 

J 0 J R d te[o,T] 

This completes the proof of the local integrability of e tC r]o with respect to H(af) 
(see Definition 12.131) . Because e tc is a Co-semigroup on H(af), a direct calculation 
shows 

e tc \d t + C*)ip(t) = d t (e tc * ip(t)) for every ip(t) £ C c °°([0, 00 ) x R d ). 

This relation together with the fact that (dt + £*)ip(t) £ C c ([0,oo),H(ao)) mean 
that both e tC ( d t + C*)ip(t) and rjo(d t (e tC ip(t))) are integrable with respect to t. 
Finally, we apply the Hille theorem (see DU771 Ch. 2, Thm. 6 ]) to obtain 

J Vo(d t (e tC *ip(t)))dt = rio(J d t (e tC *ip(t)) dt) = rjo(-ip( 0 )). 

Hence e tc r]o satisfies the conditions of Definition 13.21 □ 
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Our next goal is to show that a solution to problem £]} is unique under a suitable 
assumption of continuity. 

Definition 3.5. Let 77 be a broad-sense stationary random field with spectral mea¬ 
sure a. A family At ?? of transformed broad-sense stationary random fields is called 
strongly continuous if for every <p £ H(a) the mapping t <—> Af(p is continuous. 

Remark 3.6. ft follows directly from the above definition, that for every broad-sense 
stationary random field i], if At is a Co-semigroup on H(a ), then At rj is strongly 
continuous. 


Definition 3.7. We say that ip n —»• ip in C c ([0, 00 ), iJ(ero)) if there exists a compact 
set K such that |J n supp ip n U supp ip C K and 

lim sup || ip n (t) - ip(t)\\ H (a 0 ) = 0 . 

Lemma 3.8. If ip n —> ip in C c ([0, 00 ), H(ao)) and u(t) is a strongly continuous 
family of random fields, then 

lim sup \\u(t)(ip n (t) - ip(t))\\ L 2 , n) = 0 . 
n ^ co te[o,oo) 

Proof. Denote u(t) = Atgo- Since 1 1 —> Af.p> is continuous for every p, it is bounded 
on compact sets. Thus, from the Banach-Steinhaus theorem applied to the family 
At , we obtain that for every compact set K 

sup ||A*||iT ( <j o) —*-ir ( <To) = sup{||^V|| H(ao ) : |MU(a 0 ) = 1} = C K < 00 . 

teK 

Therefore 


lim sup \\u(t)(ip n (t) - if{t))\\ L 2 {n ) 

0 ,oo) 

< lim sup \\A*(ip n (t) ~ 1p(t))\\ H (a 0 ) 
n_>0 ° te[o,oo) 

< C K lim sup || ip n (t) - ip(t)\\ H {a 0 ) = 0 , 

n— ^°° tCzK 

where K D |J ra supp ip n U supp 'ip. □ 


Lemma 3.9. Let u\(t) and u 2 (t) be strongly continuous families of random fields. 
If for every vector (if 1 ,..., ip k ) £ C c ([0, oo),H(ao)) k we have 

p OO P OO 

( 6 ) / u 1 (t)(if 1 (t),...,i(> k (t))dt = / u 2 (t)(if 1 (t),...,ip k (t))dt, 

Jo Jo 

then u\(t) = u 2 (t). 


Proof. First we fix a vector (ip 1 ,..., ip k ) £ D k and a point s £ [ 0 , 00 ). For every 
j = 1 ,..., k, let us choose a sequence ip 3 nl (t) £ C c ([0, oo),H(<t 0 )) such that 


ton ii., 1 O) 

n—> 00 ’ 


p 3 , t £ [s, S + h] 
0 


Define ip 3 n h (t) = ip 3 n i(j^(t — s) + s). By the Lebesgue dominated convergence 
theorem we then obtain for i = 1,2 

p 00 ps-\-h 

lim / u i (t)(ip 1 h (t),...,ip k h (t))dt = / u l (s)(p 1 ,...,p k )dt. 

n ^-°°7 0 Js 

Therefore by ([ 6 ]) we also have 

1 r s + h , 1 r s +h 

- J u 1 (t)(p 1 , ...,p k )dt = - J u 2 (s)(<p 1 , ...,p k ) dt. 
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By strong continuity of u± and u 2 , when we pass to the limit with h —> 0 on both 
sides of the above expression, we obtain 

ui(s)(^ 1 ,. ..,</) = u 2 {s){ip 1 ,..., tp k ), 

which ends the proof, since pr ,..., p k and s were arbitrary. □ 


Theorem 3.10 (Uniqueness of solutions). Assume 0r) and 0d). If u± and u 2 
are two strongly continuous solutions to problem 0 , then u\(t) = u 2 (t) for every 


t > 0 . 


Proof. Since the space iL(ero) is separable, by Theorem 16.21 in Appendix, the image 
( d t + £*)(C)? o ([0, oo) x R d )) is dense in (7 C ([0, oo), H(ao)). Thus for every if G 
C” c ([0, cx)), H(a o)) we may find a sequence 

ifn = (dt + £*)</?„, for some tp n G C^QO.oo) x K d ), 

such that if n —> if in G c ([0, oo), H(ao))- 
Then by Lemma 13.81 

hm sup \\ui(t)(if n (t) - if(t))\\ L 2 f n) = 0 , i = 1,2. 
n ^°°te{ o,oo) 

Thus 

pOO 

(7) hm / \\ui(t)(if n {t)-if(t))\\ L 2 ( a) dt 

n—>°o 

< A (K) lim sup || Ui(t)(if n (t) - if(t)) || L 2 (n) =0, i = 1,2, 

n ^°°te[o,oo) 

where K D |J n suppif n U suppf/'. 

Let us fix a vector (if 1 , ..., if k ) G C c ([0, oo), H(ao)) k . For each of its elements 
we find a sequence C CU°([0, oo) x R d ) such that ( d t + C*)p l n —> if 1 in 

Cc([0, oo), H((Jq)). Because u\ and u 2 are both solutions to problem (0j), for every 
n we have 


pOO 

/ m (t)((d t + C*)ipl(t), ...,(d t + £*)ip k (t)) dt 
Jo 

poo 

= / u 2 (t)((d t +C*)p 1 n (t),...,(d t +C*)p k n (t))dt. 

Jo 

Using 0 we may pass to the limit on both sides of this equality in order to obtain 

poo poo 

/ U 1 (t)(ip 1 (t),...,if k (t))dt= u 2 (t)(if 1 (t),...,if k (t))dt. 

J 0 Jo 


Then it follows from Lemma T3.91 that ui(s) 


d 


u 2 (s). 


□ 


Let us conclude this section by summarizing the above results. 

Remark 3.11. Assume 01 ) and @ 3 ). Let a strongly continuous family of trans¬ 
formed broad-sense stationary random fields u(t) be a solution to problem 0 in 
the sense of Definition 13.21 Then Theorems 13.41 and 13.101 as well as Remark 13.61 
imply 

u(t) = e tc rj 0 , 

i.e. the semigroup solution e tc go is the unique strongly continuous solution of prob¬ 
lem 0 , up to the equivalence in finite-dimensional distributions. 




10 


4. Scaling limits of semigroup solutions 


In this section we study certain scaling properties of the semigroup solution to 
Cauchy problem ©, which is unique under conditions outlined in Remark Id. Ill 
To simplify the notation we introduce a family of scaling transformations v r 
defined for every r > 0 by 

v r <p(x) = r d p{rx) for every ip £ L 2 (dx), 

v r Z(A) = Z{rA) for every orthogonal random measure Z. 

We use the same symbol for different scalings of functions and measures. This is 
done in order to obtain the following identity 

( 8 ) [ T{v r p)dZ = I T{p\d(v r Z), 

J R d J R d 

which is an immediate consequence of change of variables in the definition of the 
Fourier transform. For a generalised random field £ we write fu r to mean £i / r {ip) = 
f{v r p) for every tp £ D. 

Theorem 4.1 (Scaling limit). Let Z$ be the orthogonal random measure associated 
with the initial condition r )o to problem © and do be its spectral measure. Suppose 
there are constants a, 7 £ R and an orthogonal random measure Z with a reference 
measure a such that 

(9) w-lim T^pt-Zq^Z. 

T —too 

Suppose further, there is a constant ft £ R and an operator V satisfying ©>) and 
such that H(a) is also V*-admissible, and 

(10) lim T 21 [ \T{e tTC {v T pp) — V 5 )}! 2 da 0 = 0 for every ip £ D. 

T->oo J Rd 

Then for the broad-sense stationary random field r]((p) = f Rd dZ we have 

w-lim T 1 (e rtc r\$)v T fj = e tv g. 

T^-oo 

Proof. For every p £ D we have the representation (see Theorem 12.91) 

e TtC r]o(v T f}p) = [ T{e tTC (v T f>p)} dZ 0 . 

J R d 

We rewrite this expression as the following sum 


[ F{e tTC \v Tf ,p))dZ Q = 

JR d 

= [ F{e tTC (v T pip) - v T <*{e tv ip)}dZ 0 + [ T{v T ‘*{e tv ip)} dZ 0 . 

JR d JR d 

For the first term on the right-hand side we use assumption m in order to obtain 


lim E 

T —too 


T 7 / F{e tTC ’ ( v T , ip) - v T « ( e tv * p )} dZ 0 


lR d 


= lim T 27 

T —too 


[ | T{e tTC * (v T cnp) - v T a (e tv * p)}\ 2 da 0 = 0. 
jR d 


For the second term, because of identity ©, we have 

[ T{v T «{e tV *p)}dZ Q = [ F{e tr ‘p}d(u T cZ 0 ). 

J R d J R d 
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Let us fix a vector {ip 1 ,..., ip k ) £ D k . It follows from assumption © that 

w-lim T 1 [ (_F{e"V}> ■ • • ,^{e tP V'}) d(v T *Z 0 ) 

T^-oo j Rd 

= f {T{e^p^},...,T{e^p k })dZ = e^r,{ip\...,ip k ) 

J R d 

and finally 

w-limT ,7 (e Tt ' c ?7o)z'Tf> = e tV V- □ 

T—y oo 

Remark 4.2. Notice that if for a suitable choice of a, /3, 7 either 

m'y r/ d ry 

T’vt* Z 0 = Zq 


or 

e tTC (y T mp) = VT<*(e tc *p) f° r every ip £ D, 

then one of the hypotheses © and GSD in Theorem ITTl is trivially satisfied. 

Remark 4.3. Observe that in the context of Theorem l4. li the limit family of random 
fields e tv ri is a solution to the following problem 

( d t u = Tu on [0, 00 ) x R d , 

\ «( 0 ) = 77 on R d . 


5. Examples 


5.1. Heat equation with Gaussian noise. First, we illustrate and explain our 
results in the simplest case of the heat equation 

( d t u = A u on [0,oo) x R, 


( 11 ) 


u{ 0) = 770 


on R. 


Here, A = <9 2 denotes the Laplace operator, or more precisely, its closure in L 2 {dx) 
with the domain Dom(A) = fF 2,2 (R) (the usual Sobolev space). We introduce the 
orthogonal random measure Zq associated with 770 by the means of the formula (see 
Theorem 12.91) 

rjo(p) = / Zo(dx) for all ip £ D. 

J R 

Theorem 5.1. Keeping the above notation, assume there exists a function / > 0 
such that 


( 12 ) Z 0 (A)= f ^J^)W(dx), 

J A 

where W is the white-noise orthogonal random measure (see Example I A 151) . Then 
the family of transformed broad-sense stationary random fields e tA ? 7 o is the unique 
strongly continuous solution to problem CD. up to the equivalence in finite-dimen¬ 
sional distributions. 


Proof. First we notice that A* = A, D C Dom(A) and e tA is a Co-semigroup of 
contractions on L 2 (dx). Consider the space H(ao) = H(f(x)dx). Observe that 
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(14) 


e tx2 < 1 for t > 0 , and therefore 

ll e ‘ A ^llff( ffo) = / |^{e tA V 2 }| 2 /(:r)<£r + f \ e tA p\ 2 dx 
Jr Jr 

(13) = [ \e~ tx2 \ 2 \T{ip}\ 2 f(x)dx+ J- f |e“ te 2 | 2 ||J r {^}| 2 da; 

Jr Z7r Jr 

< [ \J 7 {<p}\ 2 f{x)dx+ [ \ip\ 2 dx < |M|^ ( ). 

Jr Jr 

Thus H(ao) is e* A -invariant subspace of L 2 (dx). A similar calculation leads us to 

l|e* A ¥>- <p\\h(* 0 ) = [ \T{e tA p-p}\ 2 f{x)dx + [ \e tA p - ip\ 2 dx 
Jr Jr 

= f \e~ tx2 ~ l| 2 |^ r {v 5 } \ 2 f( x ) dx + [ \e tA ip - ip\ 2 dx. 

Jr Jr 

Using the Lebesgue dominated convergence theorem we obtain 

lim [ \e~ tx — l\ 2 \T{tp}\ 2 f(x)dx + [ |e tA <p — p\ 2 dx = 0 , 
t^o+J r Jr 

which shows that e tA is a Co-semigroup on H(ao). Therefore is A-admissible 

(cf. &>)). 

Then it follows from Remark 13.111 that the family e tA r/o is the unique strongly 
continuous solution to problem CD- up to the equivalence in finite-dimensional 
distributions. □ 

Theorem 5.2. Keeping the above notation, assume additionally 
lim T~ k/,2 f(x/VT)\x\ k — 1 for some k £ [ 0 , 1 ), 

, „ T—>oo 

(15) 

T- k/2 f(x/Vf)< w for every x and T and some constant c £ R 

and suppose u is a strongly continuous solution to problem CD- Then for the 
rescaled family of random fields 

u T (t,ip) = Tu(Tt)(v T -i/ 2 tp) 

and the family 


Jr \x\ ' 


we have w-liimr-j-oo u T (t) = w(t). 

Proof. By Theorem 15. II we have u(t) = e tC r)Q. Let us observe that for every ip £ D 

= i “p (- ffr ) 

\x - VTy\ 2 


tT A/ 


\JAir tT 


exp 


nx — y l j/| \ 
V 4 Tt J 


^X^ exp ( 


4 Tt 
\x/VT - y\ 2 
At 


v(y) dy 

'jp(y)dy = v T -i/ 2 (e tA <p)(x). 
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Thus e £A 77o trivially satisfies assumption m of Theorem ixn with a = (3 = — ^ 
(see Remark 14.21) . We also have 


rVz 0 (r l / a A) = f i T -i / 2 A (x)JJ(x)W(dx) 

JR 

= t^ J ^ f(x/Vr)w(dx/Vr) = j T- k ^\Jf(x/Vr)w(dx), 

where we used the scaling property of the white noise W (rA) = yfrW (A). Denote 
Z(A) = f A \x \~ k / 2 W(dx). Then 


E\T k r L Z 0 (T- 1 / 2 A)-Z(A)\ 2 =E [ T~ k ^f{x/Vr) - T J^ W{dx)' 2 

J A 


\ x \k/2 


= ±:\T-^f { x/VT)\x\ k - 1| 2 dx. 


Hence, using assumptions m and by the Lebesgue dominated convergence theo¬ 
rem, we obtain 

lim ElT 1 ^ 1 Z 0 (T~ 1/2 A) - Z(A)\ 2 = 0, 

T—>-oo 


which because of Proposition 12.61 leads us to 


w-limT 4 v t -i/ 2 Zq = Z. 

T —>oo 


Finally, in view of Theorem 14.II we have 
(16) w-lim u T (t,-)± [ 6 tX ^j ' U v(dx). 

T-s-oo J R \x \ k / 2 


□ 


Remark 5.3. The family obtained as the limit in the expression m is a solution 
to the problem 

( d t u = A u on [0, oo) x R, 


(17) 


*(0) = Co 


on R, 


where £o(<p) = / R \ x \\/2 W{dx). It is worth noting that if the function / in 

formula m is integrable, then the initial condition ?yo may be represented as a 
regular random field (see Subsection 12.61) . Meanwhile, since \x\~ k is not integrable, 
the same is never true for the initial condition £o of the limit problem (HZD- However, 
for t > 0, the semigroup solution of (ED may again be represented as a regular 
random field, since \x\~ k e~ 2tx is always integrable for any k g [0,1). 


5.2. Pseudo-differential operator with self-similar noise. In this example we 
consider the following problem 

( d t u = Vu on [0, oo) x R, 


(18) 

' u(0) = Co on 

Here the initial condition has the form 


k G [0,1), 


where W is the Gaussian white-noise. The operator V is pseudo-differential, namely 
E{Vv) = p(x)T{<p} 

for a given continuous function p : R —> C such that Rep < 0. 
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Theorem 5.4. The family of random fields e tv fo is the unique strongly continuous 
solution to problem (USD, up to the equivalence in finite-dimensional distributions. 

Proof. The proof is entirely analogous to the proof of Theorem 15. II 

The spectral measure of £o is \x\~ k dx. Consider the space H{\x\~ k dx). Observe 
that for t > 0 we have \e tp ^\ < 1 and then by a similar calculation as in m 
and (fT-fll we obtain that H(\x\ k dx) is (P*-admissible (cf. dp))). 

Then it follows from Remark 13.111 that the family e tv fo is the unique strongly 
continuous solution to problem (USD, up to the equivalence in finite-dimensional 
distributions. □ 


Theorem 5.5. Assume there exists a continuous real function q : R — > R. and a 
parameter a such that 

(19) lim Tp(T a x) = q(x) for every i£l 

T->oo 

and suppose u is a strongly continuous solution to problem (USD. Then for the 
rescaled family of random fields 

u T (t,ip) =T ! i i ' k ~ 1 ^u(Tt)(vT'*T) 

and the family 

ro(t)= /£d^M w 

Jr \x\ ' 

we have w-liniT-xjo u T (t) = w(t). 

Proof. Denote Z${dx) = a .|i / 2 W(dx) and cto = jy^dx. Then 

i ' t ° z ° {a) 1 L 1t ‘ a{x) w" W(ix) i L ia{x) w^ w(T ° dx) 

Hence 

T^ k -^u Ta Z 0 (A) = Z 0 {A) 

and Z 0 satisfies assumption © of Theorem 14.11 with 7 = %(k — 1) (see Remark [OP . 
Define an operator Q by J-{Qy} = q(x)p{y}. Then we have 


[ | P{e tTV (v T «y) ~ ur»(e tQ y)}\ 2 da 0 
JR 


= / e 


tTp(x) 


H V'-pOL </?} — x ^J-{i'T a y}\ 2 dtTo 


— / |giTp(T“x) _ e tq{x) I 2 


_ y-27 


'J-'Q 


\T a x\ k 


dx 


0 tTp(T a x) _ „tg(x)|2 


I^MI 


■ dx 


Therefore, because of USD and since both Rep and q are non-positive, we may use 
Lebesgue dominated convergence theorem in order to obtain 

lim T 2j [ \p{e tTV — VT«{e tQ <p )}| 2 da Q = 0. 

Jr 


This shows that assumption m is satisfied and it concludes the proof. 


□ 
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Remark 5.6. The family w(t) obtained as the limit in the previous theorem is a 
solution to the problem 

( d t u = Qu on [0, oo) x R, 

\ u( 0) = Co on R. 

Moreover, it follows from assumption TO that limT->oo Tp(T a x ) = q(x) and there¬ 
fore q(x) = q( l)|ai| —1 /“. Thus Q is in fact a fractional Laplace operator, namely 
for s = and c = —q( 1) we have Q = —c(—A) s . 


5.3. Heat equation with non-gaussian noise. In the last example we would 
like to present results similar to those obtained in [ LW981 . Consider the following 
problem 


( 20 ) 


d t u = A u on [0, oo) x R, 

u(0) = G(t 7 o) on R, 


where 770 is a regular broad-sense stationary Gaussian random field with the or¬ 
thogonal random measure Z 0 (A) = f A yjf{x) W(dx), and the function G is such 
that EG(rio(x)) 2 is finite. Then the function / is integrable, and the field G{rjo) 
may be represented in the basis of Hermite polynomials (see |Nua95l . Thm. 1.1.2]) 

OO 

G(vo{x)) = Y h n{vo(x)), 

z ' n! 

n =0 

where 

h n (vo{x)) = /•••/' e %x{ - yi+ - +Vn V/(j/i) • • -fiVn) W (dyi) ■ ■ -W ( dy n ) 

J R JR 

and 

<2i) 

(see |Nua95l Prop. 1.1.4]). This allows us to extend G(rjo) to a generalised random 
field by the following expression (cf. Subsection 12.611 


G(vo)(<fi) = Y ~1 / :F {tp}{y 1 + ...+y n )<Y) V Khz) W{dyi) 

n—0 U ' • dRn »=1 

for every ip £ D. 


Proposition 5.7. Each random field h n (rjo(x)) is broad-sense stationary with an 
orthogonal random measure 

r. n 

(22) Z n (A) = / l j4 (2/i + ... + 2/n)(g)v / 7(yObF(d 2 / l ). 

J * n i =1 

Proof. We proceed by a direct calculation. Denote y = yi + ■ ■ ■ + y n - 


Eh n (r]o(x 1 + h))h n (i] 0 (x 2 + h)) 

a n r n 

e dx 1+ h)y 0 W{dyi) / vTfo) W{d yi j) 

l n -i JR n 7 


2=1 




= / e ih "e fa i v e- ih ye- ix -y^f{y i )dy i = Eh n {r ]0 {x 1 ))h n {r ]0 {x 2 )). □ 

2=1 
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Remark 5.8. Formula (B^l) determines the spectral measure a n of the random field 
hn{r]o(x)), namely 

(r n {A) = [ 1 a(vi +■■■ + y n )f(yi)---f(y n )dy 1 ---dy n = [ f* n (y)dy. 

J R" J A 

Now we are in a position to identify the solution to problem (l20l) . 

Theorem 5.9. The family of random fields e tA G(r]o)(ip) = G(r]o)(e tA ip) is the 
unique strongly continuous solution to problem m, up to the equivalence in finite¬ 
dimensional distributions. 

Proof. Because EG(rjo(x)) 2 is finite, so is the spectral measure of the initial condi¬ 
tion. Then we proceed exactly as in the proof of Theorem 15.11 □ 


Theorem 5.10. Let u(t) = e tA G(p o) be a solution to (l^Ul) . Denote i?(|x|) = 
Erio(0)rio(x) and m = min{n : c n ^ 0}, where c n are given by (EH)- Suppose the 
following 

(1) R : [0, oo) —> [0,1] 

(2) R is non-increasing. 

(3) lim^^oo i?(a;)|a;| 1_fc = 1 for some k e [1 - 1). 

Then for the rescaled family of random fields 

u T (t,ip) = T^ ( d~ k ' ) u{Tt){v T - l/a y>) 


and the family 


>(».*» =cr' 2 ^ t / 

III'. J^m 


E{ip^e ~ t( - Vl+ - +Vm ^ 2 

\yi---y m \ k/2 


0 W(d yi ) 
2=1 


we have w-limT->oo u T (t) 


= w(t). 


Proof. Since we already established the scaling properties of the Laplace operator 
in Theorem IQ we need only to study the behaviour of the rescaled orthogonal 
random measure. 

Denote 7 = ^(1 — k). Similarly to a calculation in the proof of Theorem 15.21 we 
have 

T 7 Z n (T~ 1 ^ 2 A) = T 7-n / 4 [ 1 A (yi + ■ ■ • + Vn) 0 JfiVi/Vf) W(d Vl ). 

i= 1 

Since 

R(\x\) = [ e“ s /(£)d£ = 2 [ cos(a;£)/(£) 

JR JO 

then because of a Tauberian theorem (see [BGT891 Thm. 4.10.3]) we have 
lini f{x)\x\ k = C kt 

x —>0 

for a constant C k (see (BdT89l or ILW98] for details). Thus 

lim T~ k ^ 2 f(x/VT)\x\ k = C k and lim T~ k / 2 f(x/Vf) = 

x->0 T-¥o o \X\ k 

Hence for n = m we obtain 

r m 2 

lim E(T>Z m (T-V 2 A)-C£ /2 / l A (vi+. ■ -+y m ) 0 \yi\~ k/2 W(d yi )) = 0. 

T-J-OO V J R m ' 
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We also have /*”(£) = 2 f^° cos(fx)R(x) n dx and therefore by using Lemma HOI 
from Appendix we know that there exists S such that for every £ G [0, <5) and 
every n 

r (n+1) (0 <r n (o 

Therefore if A C (—6, 6) and for n> m 

E\Z n+1 (A)\ 2 = [ f< n+1 \y)dy< [ f* n (y) dy = E\Z n (A)\ 2 
J A J A 

or, for every bounded A and sufficiently large T, 

E\T^Z n+x {J~ Y l 2 A )\ 2 < E\T^ZniT- 1 / 2 A )\ 2 
(ra + l)! — n\ 

This gives us the following estimate 


e\ Y —T^Z n (T- xl2 A) 

I n\ 

n>m 


E -^E\T^Z n (T-V 2 A )\ 2 


< E\T^Z m+1 (T-^ 2 A)\ 2 ^ c| 

^ n\ 

n>m 


(m + 1 )! 


The series ^ 


lim 

T—>oo 


n>m HF c n convergent because of the Bessel inequality. Since 

E\T^Z m+l (T~ l / 2 A)\ 2 

(m + 1)! 


we obtain 

lim e\ y ^-T^Z n {T~ 1/2 

T—»oo I Z-^ n \ v 

n>m 



Finally, because of Proposition 12.61 we have 


w-limT 1 / 4 u(Tt)(^ T -i/2<p) 

T—>oo 


d / 2 

k m\ 


Ey}(y 1 + ...+y m )e- t ^+-+y^ 
\yi---y m \ k/2 


\W(dyi). 


i= 1 


□ 


6. Appendix 

Here we present auxiliary results required in the proof of Theorem 13.101 and a 
simple lemma we need to estimate an expression in the proof of Theorem 15.101 

Lemma 6.1. Let X be a locally convex separable metric linear space and K be 
a compact set. If Y is a dense linear subspace of X, then C(K,Y) is dense in 
C(K,X) (both with the uniform convergence topology). 

Proof. It follows from [MilOll Thm. 4.2.14] that (X, Y) is an ANR-pair in the sense 
of [MilOll p. 266]. Then by [MilOll Thm. 4.1.6] there exists a homotopy h on 
[0,1] x X such that h(0) = idx and h{t){x) G Y for every t > 0 and x G X (Y is 
homotopy dense in X). 

Let us take a function g G C(K,X) and put u t (k ) = h(t)(g(k)). Then we have 
uq = g and for every t > 0, u t G C(K, Y), and 

limsup ||u t - tto||x = 0. □ 

t-¥ 0 

Theorem 6.2. Assume @0 and (0b). The operator 

{dt + £*) : C c °°([0, oo) x R d ) C c ([ 0, oo), H(a 0 )) 

has a dense image. 
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Proof. Let g £ C c ([0, oo), H(ao)) and denote T 
formula 


i(t) = 


f 

Jo 


*(T-s)r 


g{s) ds, for t <T. 


Then we have 


supsuppg(s). Define u by the 


(d t + £*)u{t) = git) and u(t) £ C{[0,T],H(a o )), 

because u is a mild solution to an appropriate (reversed in time) inhomogenous 
initial value problem (see Paz831 Sec. 4.2]). We know that C f °°([0. T] x M d ) is 
dense in C([0, T], D). The space H{cro) is by definition a Hilbert completion of D 
(see ([3])), and hence contains it as a dense linear subspace. Since H{of) is separable, 
it follows from Lemma RTT1 that C([0,T], D) is dense in C{[0,T],H(cro)). 

We may therefore take a sequence 

u n £ C^([0,T\ x R d ) c C c °°([0,oo) x R d ), 
u n t u in C([0,T],H{ao)). 

Then lim„^ 00 (9 t + £*)u n = ( d t + C*)u = g and thus 

(d t +£*) : Cf°([0,oo) x R d ) -» C c ([0,oo),7L(a 0 )) 
has a dense image. □ 


Lemma 6.3. If R : [0, oo) —> [0,1] is a non-increasing and integrable function, 
then there exists 5 > 0 such that for every n and every £ £ [0, 5) we have 

POO POO 

/ cos{^x)R n+1 {x) dx < / cos i^x)R n {x)dx 
Jo Jo 

Proof. Since the function 1 1 —> /(* R(x)( 1 — R(x)) dx is positive and non-decreasing 
there exists such T that 

0<— f R{x){l — R{x)) dx — f R{x){l — R(x)) dx. 

2 Jo Jt 

Hence for S = ^ and every ^ £ [0, <5) 

pT pOO 

0< / cos(^a;)i?(x)(l — R(x)) dx — / | cos(^x)|i?(x)(l — R(x)). 

Jo Jt 

Now suppose that for some n and and every f £ [0, S) we have 

pT pOO 

0< / cosi^x)R{x) n {l — R{x)) dx — / | cos{l;x)\R{x) n {1 — R(x)). 

Jo Jt 

Then since R is non-increasing 

[ cos i£x)R n+1 (1 -R)dx>R{T) [ cos (£x)R n (1-R)dx 
Jo Jo 

p OO poo 

> R{T) / | cos(Cx)|J?”(1 - R)dx> / \ cos{fx)\R n+1 {l - R)dx. 

Jt Jt 

Thus by the induction principle we prove that for every n and every £ G [0, S) 

pT pOO 

0< / cos(£r)i?(x) n (l — R(x)) dx — / | cos(£x)\R(x) n (l — R(x)) 

Jo Jt 

pOO 

< / cos{fx)R{x) n {l — R(x)) dx 

Jo 

and finally 

POO POO 

/ cos(£x)R(x) n + 1 dx< cos (£x)R(x) n dx. □ 

Jo Jo 
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